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We develop a novel self-consistent approach for studying the angle resolved photoemission spec- 
tra (ARPES) of a hole in the t-J-Holstein model giving perfect agreement with numerically exact 
Diagrammatic Monte Carlo data at zero temperature for all regimes of electron-phonon coupling. 
Generalizing the approach to finite temperatures we find that the anomalous temperature depen- 
dence of the ARPES in undoped cuprates is explained by cooperative interplay of coupling of the 
hole to magnetic fluctuations and strong electron-phonon interaction. 
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Polaron is one of the fundamental problems extensively 
studied both theoretically and experimentally for a long 
term [l|. However, most of the theoretical works are still 
restricted to some limiting cases, or rely on unjustified 
approximations. This situation has been recently im- 
proved dramatically by the combination of the rapid ad- 
vances in the angle resolved photoemission spectroscopy 
and the numerically exact Diagrammatic Monte Carlo 
(DMC) simulation combined with the stochastic analytic 
continuation 

Especially the ARPES in the parent compounds of 
the high temperature superconductors, i.e., undoped 
cuprates, turned out to be an ideal arena for studying 
the dynamical properties of polaron formation of the sin- 
gle hole in the antiferromagnetic background. This is ex- 
pected from the ionic nature of the parent compounds 
and the associated strong electron-phonon interaction 
(EPI), together with the strong electron correlation as 
evidenced by the realized Mott insulating state. There- 
fore, the interplay between the magnetism and EPI is a 
key to resolve the quantum dynamics of the doped carrier 
into the cuprates. 

This problem has been theoretically studied [1, Q in 
terms of a hole in the t-J model coupled by short range 
Holstein interaction to optical phonons (t-J-Holstein 
model). The theoretical predictions were verified exper- 
imentally 0, Q and it was shown that the real quasi- 
particle peak has only a tiny weight at lower binding 
energy compared with the involving multi-phonon ex- 
citations Franck-Condon peak. The former one can be 
hardly observed because of tiny spectral weight while the 
latter broad one reproduces the dispersion of the pure t-J 
model. 

It has been shown [^, 0] that the polaronic effect is 
enhanced by the entanglement of the interactions of a 
hole with magnons and phonons, and this interplay is 
the unique feature of the cuprates. Indeed, the EPI alone 
is absolutely unable to explain the temperature depen- 
dence of ARPES because experimentally found temper- 



ature dependence of ARPES is considerably larger than 
that predicted by polaronic theory Q. A magnetic sub- 
system alone is also not a suitable candidate since the 
typical energy scale of magnons ~ 2J « 0.2eV is even 
larger than that of phonons ~ ~ 0.04eT^. Given such 
desperate situation, there is a temptation to explain the 
temperature driven peak broadening by the approach- 
ing of the system to the Neel temperature where quan- 
tum/thermal fluctuations destroy the antiferromagnetic 
background of the t-J model. Recent studies revealed 
one more puzzle of the temperature dependence ques- 
tioning the polaronic scenario 0, 0, 01- The tempera- 
ture dependence of the line-width is linear in the range 
400K< T <200K @ and extrapolates to zero line-width 
at zero temperature. From the theoretical point of view, 
it is a challenge to study the temperature dependence of 
the Lehman spectral function (LSF) for the t-J-Holstein 
model in the intermediate or strong coupling regime in a 
reliable way, which has never been achieved to the best 
of our knowledge. 

In the present Letter we solve t-J-Holstein model by a 
novel Hybrid Dynamical Momentum Average (HDMA) 
self-consistent method uniting the advantages of Momen- 
tum Average (MA) approach keeping the essential 
information on the magnon dispersion, and Dynamical 
Mean Field (DMF) technique, properly taking into ac- 
count strong but essentially local coupling to the lattice. 
Comparing results of HDMA method with exact data ob- 
tained by DMC approach we show that HDMA method 
provides accurate results for t-J-Holstein model where 
quasi-particle weakly interacts with delocalized magnons 
and is strongly coupled to local vibrations. Making a 
generalization of HDMA technique to finite temperatures 
we show that the basic features of anomalous tempera- 
ture dependence of ARPES in undoped cuprates can be 
explained by mutual interplay of magnetic and lattice 
systems in the t-J-Holstein model. 

The Hamiltonian of the t-J-Holstein model in the spin- 
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wave approximation 3, ll|, reads 

H = ujo'Y^ blbk + gLUQ ^ ^hlhk^qbq + H.c. 
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[Mk,qhihk-qaq + H.C. (1) 



ijJkQ'k'^k T- / ^ \ ^^j^k,qi>'k"'i^—q"-q 

k k.q 

where h\.^ and 5^ are the creation operators of a 
hole, a magnon and a phonon of momentum fc, respec- 
tively. The hole motion is associated with the creation 
and annihilation of magnons of energy ojk = 2 Jy^l — 7^ 
ilk = (cosfc X + cosfcy)/2) with coupling vertex Mk.q — 

U(Uq-1k-q 

V. 



Vq"fk)/vN where Uq 



q - -sgn{'^q)yj{l - a,)/(2aq), ttq = Wq/2J, and N is 



the number of lattice sites. The short-range interaction 
between the hole and local distortions due to dispersion- 
less optical vibrations with frequency is described by 
the coupling constant g. For the following we use the cor- 
responding to experiment values J/t = 0.3, (jJ^/t = 0.1 
0, measure all energies in units of t and assume Planck 
and Boltzmann constants equal to unity. 

The generic features of the model ([T]), causing difhcul- 
ties to semi-analytic approaches, is the intrinsic interplay 
between interaction of a hole with magnons reducing the 
spectral weight of its quasi-particle as well as reducing 
its bandwidth and coupling to local phonons backing the 
self-trapping of the quasi-particle. Brutal force disentan- 
gling of these two contributions is impossible because the 
energy scales of two processes are of the same order [13 1 
and the only attempts which were successful so far, in 
quantitative description of the spectral properties of the 
model ([T]) at zero temperature, were based on numerically 
involved methods, such as exact diagonalization 1^ 3 
or DMC y techniques. However, the results of the for- 
mer method are limited to small lattices while the latter 
one, working in the thermodynamic limit, requires ex- 
tremely extensive numerics efforts at finite temperature 
due to the "ill posed" nature of the analytic continuation 

In spite of the same energy scales involved into the 
hole-magnon and hole-phonon couplings, these two in- 
teractions are profoundly different since the coupling 
to magnons is essentially momentum dependent and al- 
ways weak whereas that to phonons is local and can be 
strong. Indeed, spin S — 1/2 cannot flip more than 
one time assuring that each site can not possess more 
than one magnon [l3| and, thus, the weak-coupling Self- 
Consistent Born Approximation (SCBA) is satisfactory 
for small values of J/t 0, [i3]- To the contrary, SCBA 
fails for EPI even in the intermediate coupling limit 
Therefore, to cope with t-J-Holstein problem it is enough 
to treat the essential momentum dependent coupling to 
magnons within the SCBA and to sum vibrational vari- 
ables nonperturbatively, at least in some local approxi- 
mation. Nonperturbative local approaches, valid at any 
coupling strength and neglecting the fc-dependence in the 
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FIG. 1: The ground state energy (Egs) and the spectral 
weight {zgs) as a function of the EPI strength A = g'^ujo/4:t in 
the DMC approach (line) and in the HDMA method (points). 



self-energy E/i_p/j(fc, w), are DMF technique [l8| and re- 
cently developed MA method [9] providing explicit form 
for the hole self-energy due to hole-phonon interaction in 
terms of a continued fraction 
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The difference in DMF and MA lies in the definition of 
a{u!) which is a function that has to be fixed by a self- 
consistent procedure in the DMF approach while is iden- 
tified with the fe-average of the bare Green's function 
in the MA scheme. Obviously, MA scheme is preferable 
when one is interested in properties of 2D model ([T]) with 
highly anisotropic coupling Mk.q. 

Summarizing the above considerations a reasonable 
self-consistent procedure expresses the total self-energy 
of the hole as the the sum of the self-energies caused by 
magnetic and phonon subsystems 

i:tjHik,u;) = Ef_<^^^(fc,c^) + Eh_ph[«tj^^(c^)]. (3) 

Weak and highly anisotropic interaction with magnons is 
taken into account in the SCBA 
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and the Q!(a;)-function for hole-phonon self-energy 
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is expressed in terms of momentum average of "bare" 
Green function whose k-dependence is determined by the 
hole-magnon self-energy ([4]) in the SCBA. The equations 
([2][5|) constitute the self-consistent set of the HDMA ap- 
proach and can be solved by standard iterative methods 
on a finite lattice by evaluating a; on a finite mesh of 
points. 

The set of eqs. (I2][5]) has the same structure as ob- 
tained in the DMF formulation of the t- J-Holstein model 
[l9| . with the important exception that the a(a;)-function 
is determined not from the purely local self-consistent 
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FIG. 2: The LSF A(ti/2,ti/2) for different EPI A in the 
HDMA (solid hne) and DMC (dotted line) approaches. Ver- 
tical line in the panel with A = 0.462 indicates the position of 
the ground state quasi-particle peak. 



DMF condition but defined through the momentum av- 
erage of the "bare" Green function containing the 
anisotropic self-energy of two dimensional t-J modeL 
Within the framework of DMF approach the t-J-Holstein 
model is indistinguishable from t-J, model where the 
hole coherent motion is suppressed [l9|. To the contrary, 
HDMA approach preserves coherent motion of the hole. 
The ground state energy, Eqs, and its spectral weight, 
ZGS = (l - dT,tjHldu!\^^^^^) are in good agreement 
with the data of numerically exact DMC approach Q 
and the crossover to the strong coupUng Hmit at A > 0.4 
is perfectly reproduced (see Fig. 1). 

For further check of validity of the HDMA scheme, 
we compare the spectral function calculated on a lat- 
tice 64 X 64 with that obtained by the approximation- 
free DMC technique for pure t-J model (A = 0) and t- 
J-Holstein model in the weak (A — 0.1), intermediate 
(A = 0.289), and strong (A = 0.462) couphng regimes 
(Fig. 2). The very good agreement of the overall shapes 
is observed for all coupling regimes and mismatch of fine 
details can be attributed to finite size effects of 64 x 64 
lattice [27II and the "local" approximation used in the 
present approach [eq. ^] that, at strong coupling, gives 
the typical oscillations with the period of phonon energy. 
Apart from these details our approach is reliable in all 
coupling regimes and gives the spectral function with a 
computational effort much less than spent by DMC. 

An important advantage of our scheme is that its gen- 
eralization to finite temperature is straightforward. Per- 
forming analytical continuation of S^_'^^^(A;, w) to Mat- 



subara formalism one gets (20|, [21 1 



n,f,(Wq)) 



E 



(6) 



where nb(w) is the Bose- Einstein factor. For the general- 
ization of Eq. ([2]) one notes that hole self-energy due to 
hole-phonon interaction Y,h-ph{^) for the model where 
the hole can interact with phonon only when it is on the 
site i the temperature dependence can be included in the 
exact way giving i23i] 
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where x = cxp(— /Jwo) and 



, (ti-2)(s.^o)^q(c^ + 2i>,o)q(^ + 3J^ 
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Since the a(a;)-fimction in Eq. ([5]) is reduced to a local 
momentum independent value, above expression is also 
valid for our scheme. 
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FIG. 3: The LSF (A(7r/2, 7r/2)) for different EPI couplings A 
and different temperatures /? = t/T. 
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The Eqs. (|3I5I6I7P provide a set of self-consistent equa- 
tions that can be solved iteratively typically within 40 
iterations. We verified the relevance of the relation ([7]) 
for our scheme checking the sum rules [§] and found that 
the first three sum rules for the LSF are satisfied at any 
temperature and EPI with high accuracy. Therefore, our 
results for peak energy and linewidth, determined mainly 
by the first and second sum rule, do not lose accuracy 
from the approximations made to obtain the HDMA self- 
consistent scheme. 

The temperature dependence of LSF at fc = (7r/2, 7r/2) 
at different values of EPI is shown in Fig. 3. The trends 
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FIG. 4: The peak's width, Alu, as a function of the temper- 
ature, T, for the t-J model (A = 0) and in strong EPI limit 
(A = 0.462 and A = 1). The temperature T is defined in units 
of Kelvin assuming t = 0.4eV. 



in the t-J Holstein model we have shown that the ori- 
gin of the anomalously large temperature dependence of 
the ARPES in the undoped parent compound of high 
temperature superconductors originates from the con- 
structive interplay between magnetic and strong electron- 
phonon interactions. 
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for peak position and linewidth are in agreement with ex- 
perimental data [1, H, [i^, HI] . With the increase of tem- 
perature the binding energy and width of the main broad 
peak increase while its intensity decreases. It is seen that 
all temperature driven effects are more pronounced in the 
t-J-Holstein model than in the t-J model supporting the 
statement 0] that the entanglement of the magnetic 
and vibrational fluctuations is essential for cuprates and 
crucial for description of anomalously enhanced temper- 
ature driven effects in undoped compounds [8]. 

The temperature dependence of the peak width esti- 
mated through a Gaussian fitting is shown in Fig. 4. It is 
remarkable that in the strong coupling regime of the EPI 
the peak width is almost constant up to T ~ Li;o/2 ~ 
200K and then, for T ^ wo/2, demonstrates linear de- 
pendence which, in according with experiment, can be 
naively extrapolated to almost zero value at zero tem- 
perature. Note, the temperature dependence is strongly 
enhanced in the strong coupling limit A > 0.4 of the 
t-J-Holstein model giving in this limit, in contrast with 
purely polaronic or purely magnetic models, correct order 
of magnitude of the effect and even showinga good semi- 
quantitative agreement with experiment [a, @, 24, |25| . 
For A = 0.462 thepeak width, in quantitative agreement 
with experiment [61], doubles in the range from 2GGK to 
400K though the absolute value of the peak width is a 
factor of 2 below the experimental values. On the other 
hand, the absolute value of the linewidth fits the exper- 
iment for A = 1 but the enhancement of the width is a 
factor of 1.5 below that found experimentally [ij. The 
above discrepancies can be attributed to the fact that 
the longer range hoppings of more realistic tt't"- J model 
are missing in the t-J model or as well to the fact that 
in the realistic systems the holes are coupled to several 
phonon modes through the EPI of different strength [2^ . 

In conclusion, by using a new hybrid dynamic momen- 
tum average approach to the calculation of a hole LSF 
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